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Abstract 

The control of quantum phenomena is a topic that has carried out 
many challenging problems. Among others, the Hamiltonian identifica- 
tion, i.e, the inverse problem associated with the unknown features of a 
quantum system is still an open issue. In this work, we present an algo- 
rithm that enables to design a set of selective laser fields that can be used, 
in a second stage, to identify unknown parameters of quantum systems. 

1 Introduction 

The possibility to use coherent light to manipulate molecular systems at the 
nanoscale has been demonstrated both theoretically [T] and experimentally [TS] . 
Different types of methods have proven their relevancy for various settings, rang- 
ing from electron to large polyatomic molecules [51 [SI [HI [HI [Sj ■ 
At the same time, the ability to generate a large amount of quantum dynamics 
data in a small time frame can also be used to extract from experiments the 
values of unknown parameters of quantum systems. The corresponding inverse 
problem, usually called Hamiltonian identification has recently been subject to 
significant developments through encouraging experimental results [4]. 
Various formulations in an optimization settings have been studied. Because of 
the nature of the available data, zero order methods were first tested, see e.g. 
the technique of map inversion [TB] . The use of optimal control techniques was 
then introduced [31 [S] . 

Contrary to this last class of methods, we present in this work a methodology 
that enables to handle situations where the experimental measurements are pro- 
vided only at a given time. Our approach is based on a precomputation that 
provides a family of selective laser fields. Roughly speaking, these laser fields are 
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designed iteratively to highlight variations in the parameters that are subject to 
the identification. In a second stage, these fields and the experimental measure- 
ments are used to assemble a nonlinear system satisfied by the to-be-identified 
parameters. 

The paper is organized as follows: the optimization framework and the assump- 
tions we use are presented in Sec. [5] In Sec. [31 the structure of our algorithm 
is given. The procedures used in the two parts of this algorithm are described 
in Sec. 2] and Sec. O The identification step is explained in Sec. [SI Details 
about practical implementation and some numerical results are given Sec. [71 
We conclude with some remarks in Sec. [SI 

Throughout this paper, is a spacial domain in M.'^, d — 1,2, 3, denotes the 
space of complex valued square integrable functions over fl, and (., .) the usual 
Hermitian product associated to . The following standard convention is used: 



the set of all linear operator from into L^. Finally, we use 3?(z),3(z) to 
denote respectively the real and the imaginary part of a complex number z. 

2 The identification problem 

We first introduce the model and the framework used in this paper. 
2.1 Control of the Schrodinger Equation 

Consider a quantum system G H^, with norm HV'IIl^ = 1, evolving according 
to the Schrodinger equation 



where Hq is the kinetic energy operator, V G C{L^;L^) the potential operator 
and /i G C{L^;L'^) the dipole moment operator coupling the system to a time- 
dependent external laser field e{t). In this context, £ reads as a control since it 
can be chosen by the experimenter. 

In the settings we consider here, we assume that the internal Hamiltonian 
H = Hq + is known so that the goal is to identify the dipole moment operator 
fi. The generalization to the identification of V should not give rise to any 
particular problem and is left to a future contribution. 

The basic hypothesis made on /i is that it belongs to (or actually can be 
conveniently approximated by) a finite dimensional space spanned by some basis 
set Bf_, = (/i£)£=i,...,L. 



{a\0\b) := (a, 0{b)), aEL^, b E , Oe C{L^;L^) 



{ 




[Ho + V + £[t)y\i^ 



(1) 



2.2 Experimental measurements and controllability 

In order to perform the identification, we assume that given a time T and a laser 
field £ € L^(0,T), the experimenter can measure, for some fixed state & L"^, 
with norm HV'ilIt^ = 1, the value <f{n,s) := {iIji,4>{T)). 

Note that all what follows still holds when considering several measurements a 
time T, i.e., in the case where a set of measurement ((V'^, V'(^)))^=i...,p) with 
p > 1, is known. 

Finally, we assume that the system under consideration is wavefunction control- 
lable, i.e., that e G L^{0,T) i— > ip{T) is surjective. 



2.3 Formulation of problem 

Our identification method is based on a particular formulation of the identifi- 
cation problem that we now briefly introduce. 

Denote by jj,* the actual dipole moment operator of a given system. The solution 
= yu* of our problem also solves the minimization problem: 

ri^lr2, -<^(M*,e)P- (2) 

This settings highlights the fact that as long as /i 7^ /i*, a selective laser field 
should be designed so that the difference between ij, and fx* is discerned through 
the measurement </?(/L<,e). 



3 Structure of the algorithm 

Our algorithm consists in designing, through a finite iterative procedure, a set 
of selective laser fields. We start with the general structure of our algorithm. 
Details about its steps are given in the next sections. 

3.1 The selective laser fields computation greedy algo- 
rithm 

Starting from the basis set Bfj, = {i^i)i=i,...,L, the algorithm builds up iteratively 
a set of L selective laser fields as follows. 



Algorithm 1 (Selective laser fields computation greedy algorithm) Let us define 
a laser field that solves the problem: 

sup \ip{lJ,i,s)f. 

Suppose now that at the step k, with 1 < k < L, a laser fi,eld e*^^^ is given. The 
computation of is performed according to the two following sub-steps: 



1. Fitting step : Find (ct^)j^i,...,fc-i that solves the problem: 



(3) 



m f/ie minimum mean square error sense. 



2. Discriminatory step : Find that solves the problem: 



k-l 



= argmax^^^2^Q^T,)\Lp(p^ ,e) - LpC^ a] n\ e)]'^ . 



The initialization of the algorithm is somehow arbitrary, the only require- 
ment is that has a link with the type of measurement. In our case, we decide 
to maximize it. 

Remark 1 Note that, in opposition to usual approaches (see e.g. our 
method plays the role of a precomputation step since the actual measurements 
ip{pL*,e) are not required at this stage. 

3.2 Intuitive interpretation of the algorithm 

In the first sub-step of an iteration of Algorithm [1] one looks for a defect of 
selectivity of the current laser fields e^, . . . , e*"'"^: in the case the minimum 
reaches zero, two distinct dipole moment operators give rise to two identical 
measurements when exited with the laser fields e^, . . . ,e^~^. On the contrary, 
the second sub-step aims at computing a laser field that compensates this defect. 
These two sub-steps corresponds respectively to the minimization part and to 
the maximization part of the formulation 

Remark 2 Even if no hierarchy is assumed in the basis B^^, this algorithm 
should be viewed as a first step towards future works that handle infinite dimen- 
sional systems. In such a framework, the sum '^j/^"' "would read as an 
asymptotic expansion of the dipole moment operator. 

This algorithm belongs to the class of greedy algorithms, since it follows the 
problem-solving's heuristic of making the locally optimal choice (in the second 
sub-step) at each stage with the hope of finding the global optimum that solves 



4 Fitting step 



Let us first focus on the first sub-step of the algorithm. Consider an integer k 
such that 1 < k < L and denote by K'' the functional (defined on R'^^^): 

fc-1 fc-i 

rn — l J — 1 

During this sub-step, one has to find the minimum of the cost functional K''. 
To do this, a standard global minimization algorithm applied to this minimum 
mean square error associated problem. 

Note that, for small values of L, the gradient of the functional K'' can be com- 
puted thanks to the formula: 

fc-1 

E 23? (T) - Ve'™ {T), Vi) W?™ (T), V-i)) , 



where ipfm and ijj^m are the solutions of Eq. ^ with e = as laser field. 



and fi = X]j=i CKjA*"' and 11 = 11'' respectively as dipole moment operator. The 
variation Sip°' is computed thanks to: 



-[Ho + V + e'^-^t) -Ii' '5a,/i^)]^?" 



In this way the computation of the components of Vif'^ (a) can be parallelized 
to make the use of gradient methods feasible. 



5 Discriminatory step 

To achieve the second sub-step of Algorithm [1] we adapt an efficient strategy 
usually used in in quantum control. This strategy has given rise to a large class 
of algorithms often called "monotonic schemes". For a general presentation of 
these algorithms, we refer to [IT]. 

5.1 Improvement of the selectivity of a given laser field 

Let us present in more details how this strategy applies in our case. Note first 
that, given a laser field e £ L^(0, T), and two dipole moment operators Jl and 
/X, one has: 



where O^^ — ipi.tpf, tjj and •ip are the solutions of Eq. ^ with respectively 
/J, = Ji and jj, — p. as dipole moment operator. 

In order to compare the selectivity of s and e', we introduce the functional: 

J(e) = {^{T)-i^{T)\0MT)-${T))-f3 e'{t)dt, 

Jo 

which has to be maximized. For sake of simplicity, we omit the dependence of 
J with /i and fl in the notations. 

The additional term f3 J^^ e'^{t)dt, is introduced for two complementary reasons: 
first, as it penalizes the L^-norm of the laser field, it enables to obtain feasible 
laser fields and secondly, it improves the convergence of Algorithm [5] below. 
Consider now another laser field e' G i^(0,T), and denote by ip' and 'ip' the 
corresponding solutions of Eq. ([1]) with ii — Jl and fj, = ^ respectively. We 
introduce the two adjoints states defined by: 



( IX ^ [Ho + V + e{t)Jl]x 

{ X{T) = O^, (^(T) - V^(T)) 

and 

r z£ - [H^ + V + e{mx 

\ X{T) = O^, (^(T) - V7(T)) 



(4) 



(5) 



One has: 



J{e') - J{e) = (Si^'iT) ~ 5^{T)\0^, \5^'{T) - S^jiT)) 



T 

12/ 



-P / e'\t)-e\t)dt 
Jo 

= {5^/(T) - Si^{T)\0^, - S^{T)) 

+ f {e\t)~e{t)) 

JO 

(23(xWlMl?W)-23(xWlMl?(i))-/3(e'W+£(i)))rfi, 



(6) 

where we denote 5ip'{T) = ii'{T)~$'{T) and S^/j{T) = ^{T)~i}{T). Identity ^ 
gives a criterion to guarantee that e' is more selective than e. Indeed, suppose 
that e' satisfies for all t e [0,r] the condition: 

{e'{t)-e{t)) (23(x(t)|MlV^'(i)) -23(x(i)|Ai|?W) 

-/3(£'(t)+£(i))) >0, (7) 



then J(e') > J(e). 

Various ways to ensure that ([7]) holds. For example 14], one can define e' at 
each time t as the solution of the equation: 

e'(t) - e{t) = I (23(x(t)|Ml?W) - 2S(xW|Ai|^'(t)) 

(8) 

where is a given strictly positive number. In this case, one has: 
J{s')-Jie) = {H'{T)-5iAT)\O^AH'{T)~54'{T)) 
+^J^ {e'{t)-e{t)fdt>0, 

which is the desired conclusion. In Sec. 17.11 we present an alternative that can 
be obtained in a time discretized settings. 

5.2 Discriminatory sub-algorithm 

We derive form the previous considerations the following iterative procedure to 
define a laser field e*^ that maximizes J(e): 

Algorithm 2 (Discriminatory sub-algorithm) Let Tol be a positive number. 
Consider an initial guess £q and compute the corresponding solutions of Eq. ([1]) 
with Jl and /i, say ipQ and ijjQ. Set err ~ 2.Tol. 
While err > Tol, do: 

1. Use Eqs. ^^{Ej) with e — e^, ip = and i/; = il^i, to compute xi o,n-d Xtr 
respectively. 

2. Compute simultaneously the laser field ££_^_i and the states o^nd i^i+i 
the solutions of coupled system composed of Eq. ^ with x ^ Xt, X = Xi 
and Eq. ^ with fJ- = Jl and fJ- = respectively. 

In [7], one shows that Eq. ^ has a solution and presents some efficient 
numerical nonlinear solvers to compute it. 



6 Identification procedure 

Once the L selective fields , ...,e^ have been computed, one can use them ex- 
perimentally to obtain the corresponding measurements (p{fJ,* , £^), (^(m*, £^). 



The identification procedure consists then in finding the linear combination 
{a^, ...,a^) that solves the following nonlinear system: 

in the mean sc^uare sense. In this view, the standard global optimization pro- 
cedure used for the first sub-step of algorithm can be applied to the associated 
problem. 

Note that, in a finite-dimensional settings, the existence of a solution is guar- 
anteed. 



7 Numerical implementation and results 

We give here details about the practical implementation of Algorithm [U and 
show its efficiency on an example. 

7.1 Numerical solvers 

In order to solve numerically Eq. ((1]), we use the second order Strang operator 
splitting [12]. Given M > 0, a time step At such that M.At — T and an 
approximation ipj of %l}{j./S.t) with j < M, this method leads in our case to the 
following iteration: 

V'.+i = e'^^e^^^^^*e^^^V,. (10) 
In the second sub-step of Algorithm^ Discriminatory sub-algorithm[2]is adapted 
to this discrete settings. In this way, we consider the time-discretized version of 
the cost functional J: 

M-l 

where e e R^-^-^ Fix now two discrete laser fields e and e', one can then repeat 
the computation done in Sec. IS.ll to obtain: 

M-l 



J=0 



(23(xj|AiAt(4>£j)lV'j) - 23(xj|AiAt(4.£j)IV^) 

-/3(e;+e,)), (11) 



where the vectors x, £, ip' and ip' are computed using the iteration ^TU^ with 
II = Jl/^^.tiSj, £j) and fi — /iAt(£j , Ej)- These matrices are the approximations of 
ju and /i respectively defined by: 



AiAt(£',ej) 



e 



ie'-pAt AcipAt 

= zA^(s;-s,) 

For the sake of simphcity, instead of solving the discrete version of Eq. ^ , 
we compute e'j using one step of a Newton optimization method applied to 
its corresponding term in the sum of Eq. (jlip . This strategy, and the one 
corresponding to Eq. ^ are presented in more details in [7] . Their convergence 
are proven in [lOj . 

7.2 Numerical test 
7.2.1 Settings 

To illustrate the ability of our approach, we consider a simple finite dimensional 
settings where Hq, V and /i are 3x3 Hermitian matrices with entries in C and 
S C'^. The internal Hamiltonian we consider is: 



Since Eq. ([T]) with such an internal Hamiltonian is generically controllable, we 
choose to define the basis randomly so that the systems handled by our 
algorithm are almost surely controllable. 

In order to work in a general framework, we chose /z* also randomly. In our 
example, we consider: 

2.4154 1.9335 1.5822 
fi* = I 1.9335 1.4366 1.5991 




1.5822 1.5991 1.9843 



The states tpo and ipi are 

-00 





We choose T = 40007r, which corresponds to 20 periods of the transition asso- 
ciated to the smallest frequency of the system. 



7.2.2 Algorithm parameters 



The minimum mean square error problems Q-® are solved by standard pseudo- 
Newton solvers. In order to make a global search, we repeat the minimization 
10 times with random initialization. The parameter (3 is adapted to make Al- 
gorithm [5] converge. In our case (3 = 10^^. 

7.2.3 Numerical results 

The precomputation is achieved by our algorithm in approximately 80 min CPU. 
The dipole moment operator is regained with a relative error 

^*^„ « 9.896Ge - 04, 

IIm*I|2 

in approximately 10 min CPU. The selective fields that have been obtained are 
depicted in Fig. [TJ 




14000 



Figure 1: Selective laser fields obtained by Algorithm [2l 



8 Concluding remarks 

The Selective laser fields computation greedy algorithm presented in this paper 
shows a good efficiency in a general settings. However, there is some room 
for improvement of our strategy. First, the choice of the basis B^^ could be 
improved, e.g. through an iterative procedure. Secondly, the experimental 
measurements could be used during the computation of the selective fields in 
order to design an online procedure. Lastly, some work has to be done to 
design a more specific approach to treat the first sub-step of the algorithm. The 
identification procedure presented in Sec. [6] would also certainly take advantage 
of such a study. 
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